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We examined a two-dimensional Heisenberg model with two kinds of ex- 
change energies, Jg and Jc- This model describes localized spins at vanadium 
ions in a layer of CaV409, for which a spin gap is found by a recent exper- 
iment. Comparing the high temperature expansion of the magnetic suscep- 
tibility to experimental data, we determined the exchange energies as Jg — 
610 K and Jc — 150 K. By the numerical diagonalization we estimated the 
spin gap as A ~ 0.2Je ~ 120 K, which consists with the experimental value 
107 K. Frustration by finite Jc enhances the spin gap. 
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Low- dimensional antiferromagnets with spin gap attract much interest 
due to the possible relevance to the high-T^ superconductivity. For one di- 
mension, many systems with spin gap have been examined theoretically. For 
example, the ladder model opens a spin gap up to about 0.5J. |^, The 
Majumdar- Ghosh model, a typical frustrated system, has a spin gap of about 
0.24 J. 1^, H, 1^, H A series of Heisenberg models with linearly decreasing 
exchange interactions also have finite spin gaps in the range between 0.25J 
and 0.75J. |P, [l^] In the above, J is a typical exchange energy included in 



each model. Experimentally, finite spin gaps are found in (VO)2P207 [|TT 



and SrCu203. ||T2|, |T3| These materials are represented well as ladder spin 
systems arranged in parallel and coupled weakly. 

For two dimensions, Taniguchi et al. [0] recently found a finite spin gap 
for a layered material CaV409. They estimated the spin gap as A ~ 107 
K by measuring the magnetic susceptibility of d-electron spins at vanadium 
ions (V-spins) and the spin-lattice relaxation rate of ^^V nuclear moment. 
As long as we know, this is the first experiment showing clearly the spin gap 
for (quasi-) two-dimensional spin systems. The lattice structure of CaV409 
shows that there are two kinds of important exchange interactions between 
edge sharing V-spins and between corner sharing V-spins in a layer. We 
denote the corresponding exchange energies as Jg and Jc, respectively. Katoh 



and Imada ||T^ examined the spin gap by assuming = 0. They estimated 
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as A ~ 0.1 IJe by the Quantum Monte Cairo simulation and explained the 
ground state character by a perturbation calculation, although they did not 
estimate effects of Jc- Ueda et. al. [TBI carried out a similar perturbation 
calculation for finite J^. However, since these perturbation calculations are 
not justified for realistic parameter values, their results have not been definite 
yet. 

In this letter, we construct a reasonable spin model with realistic values 
for Je and Jc and explain the measured spin gap for V-spins in CaV409. 
The exchange energies are determined by carrying out the high temperature 
expansion (HTE) of the magnetic susceptibility for the model and by com- 
paring them to experimental data. We examine the spin gap by numerically 
diagonalizing the Hamiltonian for finite systems. After the extrapolation to 
an infinite system is taken, the estimated spin gap is shown to be fairly close 
to the experimentally obtained spin gap. 

To construct the Hamiltonian, we survey the structure of a layer in 
CaV409. [l^, [T^] A layer consists of VO5 pyramid-shaped clusters with api- 
cal oxygen atoms above and below a basal plane. Each pyramid contains a 
vanadium ion roughly in its center. An electron in the de^y orbital of 
a V''"'' ion forms a localized spin (V-spin). [|18| Superexchange interactions 
may occur between V-spins in edge sharing pyramids and between V-spins 
in corner sharing pyramids by hybridization of de^y orbitals with or Py 



3 



orbitals of adjacent oxygens. These superexchange interactions contribute 
to Jg and Jc respectively. There may be also a direct exchange interaction 
between V-spins in edge sharing pyramids due to the overlap of the de^y or- 
bitals, which contributes to Jg- Since it is difficult to calculate the values of 
Je and Jc by starting from the first principle, we determine them by compar- 
ing experimental data of the magnetic susceptibility to the HTE calculation, 
as will be shown. Thus we describe magnetic properties of CaV409 by a two- 
dimensional antiferromagnetic Heisenberg model with two kinds of exchange 
energies Jg and Jc- 

The Hamiltonian for CaV409 is written as 

H^Ho + Ha + Hb, (1) 

Hq — Jc X! ' Ha — Jc ^ Si ■ Sj, Hb — Jc X! ^" ' (2) 

<i,a> <i,j> <a,/3> 

where Si {So) is the V-spin at site i (a) belonging to the A {B) sublattice. 
The lattice structure in a layer is shown in Fig. 1(a). A small circle (square) 
represents V-spins above (below) the basal plane. We have called the set of 
sites denoted by small circles (squares) the A (B) sublattice. The exchange 
interactions for {i,a) is denoted by bold sohd hne in Fig. 1 (a); the corre- 
sponding exchange energy is Jc- The exchange interactions for {{a, (3)) 
are denoted by thin (dashed) solid lines; the exchange energy is Je- We note 
that the sub-Hamiltonian Ha (Hb) consists only of spins belonging to the 
A (B) sublattice, while Hq consists of spins belonging to both. In the case 
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of Jc = 0, the Hamiltonian reduces simply to Hq. The lattice for Hq given 
by thin solid lines in Fig. 1(a) is topologically the same as the lattice of 
Fig. 2. In the case of Je = 0, spins on the A sublattice are described only 
by the sub-Hamiltonian Ha and do not interact with spins on the B sublat- 
tice, which are described only by Hb- Hence we can consider Ha and Hb 
separately. The sublattice for Ha (Hb) represented by thin solid (dashed) 
lines in Fig. 1(a) is topologically the same as the lattice shown in Fig. 2. 
Thus all sub-Hamiltonians Hq, Ha and Hb are equivalently represented by 
the unfrustrated lattice of Fig. 2. The ratio 7 = Jc/ Je changes the strength 
of frustration; the system is unfrustrated both in the limits of 7 = and 

7 = CX). 

To determine the exchange energies, Je and Jc, we carry out the HTE for 
magnetic susceptibility and compare the result to experimental data. The 



experimental susceptibility xe is shown in Fig. 3 as a function of 1/T. The 
HTE is known to precisely reproduce a high-temperature part of a thermody- 
namic quantity of quantum spin systems. JTOf We derived the HTE formula 



for the magnetic susceptibility xh of the Hamiltonian (1). It is written as 
1 3(1 + 7) Je , 3(l + 67 + 72)j2n 



XH = C 



+ ), (3) 



T 4T2 16T3 

where C = n{g^BY / '^k with n being the number of vanadium ions per gram 
in CaV409, g the (7-value, /ie the Bohr magneton and k the Boltzmann con- 
stant. To precisely compare the HTE susceptibility (3) to the experimental 



data, we introduce the quantity, xh = (xh — ^)T^- Then eq. (3) gives a 
hnear function of 1/T for xh as xh = Cq + Ci/T with cq = — 1(1 + 7) Je and 
Ci = ^^(1 + 67+7^) Jg . We plotted the corresponding experimental quantities 
Xe = {xe — j^)T'^ as a function of 1/T in the inset of Fig. 3. The coefficients, 
Co and ci, are determined so that xh represents the tangential line of xe, as 
shown in the inset; i.e. cq = —2.21 emu/gK^ and ci = 693 emu/gK. From 



the coefficients, we obtain Jg ~ 610 K, — 150 K and then 7 ~ 0.25. pO 



Using these values we plotted the xh in Fig. 3 by the dashed line. The 
exchange energy of Je — 610 K seems to be fairy large in vanadium oxides 



and becomes a half of cuprates. |13| The value of Jc is smaller than Jg but 
is not negligible, so that the real system is frustrated. We should carefully 
consider the contribution of Jc when we examine the magnetic properties of 
CaV409. 

To obtain the spin gap, we numerically diagonalize the Hamiltonian (1) by 
Lanczos' method for finite systems with the periodic and/or the antiperiodic 
boundary conditions. In the case of 7 = we can use the systems with 
A^= 12, 16, 18 and 24, which are shown in Fig. 2. However in the case of 
finite 7, systems with = 16 and 24 among them only fit to the boundary 
conditions; these finite systems are shown in Fig. 1(b). This is because 
the lattice of 7 = (Fig. 2) is more symmetric than that of 7 7^ (Fig. 
1(a)). We calculate the excitation energies A from the singlet ground states 
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to the lowest triplet states both for 7 = and 7 = 0.25 and compare them. 
For the extrapolation, we assume the system-size dependence of A ~ 
The results are shown in Fig. 4. Data for 7 = fit a straight line well 
and confirms the system-size dependence. Hence the spin gap for 7 = is 
estimated as A ~ O.lSJg in the thermodynamic limit. This agrees with the 



result obtained by the Quantum Monte Carlo calculation. We expect this 
system-size dependence is correct even for small but finite 7 and apply it to 
the realistic case of 7 = 0.25. Then the spin gap for 7 = 0.25 is estimated 
as A ~ 0.2Je in the thermodynamic limit, as shown in Fig. 4. This result 
shows that frustration remarkably enhances the spin gap. Using the exchange 
energy ^610 K, the spin gap is evaluated as A ~ 120 K. This result is 
fairly close to the experimentally obtained spin gap 107 K. 

We finally discuss the origin of the spin gap. First we consider the case of 
Jc = represented by the lattice of Fig. 2. The Hamiltonian Hq consists of 
two different types of interaction bonds. An interaction bond is one of four 
bonds forming a plaquette and another connects two plaquettes. We call 
them a plaquette bond and a dimer bond, respectively. Accordingly Hq is 
written as Hq = Hop+Hoa, where i^op (-f^od) consists only of plaquette (dimer) 
bonds. We calculated the expectation values of Hop and Hq^ in the ground 
state by the numerical diagonalization. In the thermodynamic limit these 
values are (Hqp) ~ — 0.20JeA^ and (i/od) — — 0.07JeA^5 when the numerical 
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results are extrapolated as a function of Hence the portions of the 

energy gains in plaquette parts and in dimer parts are (Hqp) / (Hq) ^ 0.74 
and (Hod) / (Hq) ^ 0.26, respectively. This result agrees with the picture in 



which singlet clusters are formed in plaquette parts. |I5|, |16[ Next we consider 
the general case of Jc 7^ in terms of the spin-spin correlation function 
{Si ■ Sj). This quantity is evaluated by the numerical diagonalization and 
the extrapolation, again. At a plaquette bond we obtained (Sj ■ Sj) ~ —0.60 
for 7 = 0.25 against (Sj ■ Sj) —0.54 for 7 = 0. At a dimer bond we 
obtained (S, ■ S^) ~ -0.24 for 7 = 0.25 against (S, ■ Sj) ~ -0.39 for 7 = 0. 
This result shows that parts gaining the correlation energy move from dimer 
bonds to plaquette bonds as 7 increases. This tendency corresponds to the 
enhancement of the spin gap with frustration. 
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Figure Captions 

Fig. 1. (a) The lattice structure for V-spins in a layer of CaV409. Small 
circles (squares) represent V-spins belonging to the A {B) sublattice 
and are in pyramids above (below) the basal plane. The exchange 
energy for a bold solid line is Jg. Those for a thin solid line and a 
thin dashed line are Jc- (b) Parts of the lattice used in the numerical 
diagonalization. 

Fig. 2. The lattice representing a sub-Hamiltonian Hq, which becomes the 
total Hamiltonian when = 0. This lattice also represents each of Ha 
and Hb, which are separated from each other when Je = 0. Parts of 
the lattice used in the numerical diagonalization are also shown by the 
thin solid line. 



Fig. 3. Magnetic susceptibilities of the experiment ||14|| and the HTE calcu- 
lation. The HTE curve is determined so that xh becomes a tangential 
line of xe as shown in the inset. 

Fig. 4. System size dependence of the spin gap A for 7 = (circles) and 
0.25 (squares). The extrapolations are done with the solid and the 
dashed lines determined by the method of least squares with data of 
> 16. 
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